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ABSTRACT 

By solving radiative transfer equations, we examine three-dimensional radiative prop- 
erties of a magnetohydrodynamic accretion flow model confronting with the observed 
spectrum of Sgr A*, in the vicinity of supermassive black hole at the Galactic cen- 
tre. As a result, we find that the core of radio emission is larger than the size of 
the event horizon shadow and its peak location is shifted from the gravitational cen- 
tre. We also find that the self-absorbed synchrotron emissions by the superposition 
of thermal electrons within a few tens of the Schwartzschild radius can account for 
low-frequency spectra below the critical frequency Vc ~ 10^^ Hz. Above the critical 
frequency, the synchrotron self-Compton emission by thermal electrons can account 
for variable emissions in recent near- infrared observations. In contrast to the previous 
study by Ohsuga et al. (2005), we found that the X-ray spectra by Bremsstrahlung 
emission of thermal electrons for the different mass accretion rates can be consistent 
with both the flaring state and the quiescent state of Sgr A* observed by Chandra. 
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1 INTRODUCTION 

How the emission comes from accreting material in the 
Galactic centre (GC) is a fundamental question for un- 
derstanding the nature of mass accretion processes feeding 
a supermassive black hole (SMBH). When mass accretion 
rate is much less than the critical value, Merit = ^Edd/c^ 
where Lsdd « 1-3 x 10^** (Mbh/M©) ergs"^ is the Edding- 
ton luminosity and c is the speed of the light, the radia- 
tion loss of accreting gas is inefficient and thus most of the 
energy generated by turbulent viscosity is stored as ther- 
mal energy of the gas and is advected onto SMBH. There- 
fore, the accretion flow becomes hot and geometrically thick 
structure. This type of accretion flow is well known as an 
advection-dominated accretion flow (ADAF: Ichimaru 1977; 
Narayan & Yi 1994, 1995; Abramowicz et al. 1995) or a ra- 
diatively inefficient accretion flows (RIAF; Yuan et al 2003; 
see also Kato, Fukue, Mineshige 2008). The ADAF/RIAF 
model is quite successful in reproducing high-energy emis- 
sion of low-luminous active galactic nuclei (AGN) and the 
GC source, which is Sgr A* as a compact radio source (Bal- 
ick & Brown 1974; Bower et al. 2004; Shen et al. 2005; 
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Doeleman et al. 2008). Actually, the stellar dynamics has 
revealed the mass of SMBH in the GC, ^ 4 x 10^ M© (e.g., 
Schodel et al. 2002; Ghez et al. 2003, 2008; Gillessen et al. 
2009). It turned out that the low-luminous material at the 
GC is associated with Sgr A*, in which the luminosity is 

iSgrA* ~ 10 *I/Edd- 

After the discovery of magneto-rotational instability 
(MRI: Balbus & Hawley 1991, 1998), which can drive mag- 
netohydrodynamic (MHD) turbulence as a source of viscos- 
ity for accretion process, it seems that non-radiative MHD 
simulations of accretion flows have been accepted as a re- 
alistic model of ADAF/RIAF (Stone & Pringle 2001; Haw- 
ley & Krolik 2001). For example, many MHD studies based 
on numerical simulations have revealed a hot and geometri- 
cally thick accretion flows with a variety of complex motions 
(Matsumoto et al. 1999; Hawley 2000; Machida et al. 2000), 
outflows and jets (Igumenshchev et al. 2003; Proga & Begel- 
man 2003; Kato et al. 2004), and oscillations (Kato 2004). 
The recent detection of variability in the X-ray and near- 
infrared (NIR) emissions at Sgr A* (Baganoff et al. 2001, 
2003; Ghez et al. 2004) may be induced by such a multi- 
dimensional structure of the flow. Nonlocal nature of radia- 
tion process is essential for testing MHD models. Therefore, 
in order to clarify the structure and the time variability of 
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accretion flows, undoubtedly fuU radiation transfer treat- 
ment of MHD accretion flows in three-dimension is indis- 
pensable. 

The pioneering work for examining MHD model of ac- 
cretion flows has been done by Ohsuga, Kato, Mineshige 
(2005: hereafter OKM05). OKM05 assumed the cylindrical 
distribution of electron temperature by adopting the balance 
equation between radiative cooling and heating via Coulomb 
collision, regardless of gas temperature in MHD model. They 
reconstructed for the first time the multi-band spectrum of 
MHD model which is consistent with the observed spectra 
in flaring state of Sgr A*. However, the spectra in quiescent 
state cannot be reconstructed simultaneously in the radio 
and X-ray bands. In the context of MHD model, this is the 
issue of what makes the difference between the flaring state 
and the quiescent state in Sgr A*. We expect that determina- 
tion of the electron temperature is the key for understanding 
the occurrence of two distinct states. 

Moscibrodzka, Proga, Czerny, and Siemiginowska 
(2007: hereafter MPCS07) present the spectral feature of 
axisymmetric MHD flows by Proga & Begelman (2003). In 
contrast to OKM05, they calculate the electron temperature 
distribution by solving the heating-cooling balance equation 
at each grid point at a given time in the simulation. More- 
over, they take into account the advective energy transport 
and compressive heating in the balance equation. It turns 
out that the heating of electrons via Coulomb collision is 
not always a dominant term in the balance equation. Unfor- 
tunately, MPCS07 failed to reproduce the radio and X-ray 
spectra by emission from thermal electron. They conclude 
that a contribution of non-thermal electrons offers a much 
better representation of the spectral variability of Sgr A*. 

One difficulty in theoretical studies of radiative feature 
of hot accretion flows is that radiative energy transfer plays 
a critical role for determining the electron temperature. For 
example, radiative heating/cooling via synchrotron emis- 
sion/absorption and Compton processes may dominate com- 
pressional heating and coUisional heating in the energy bal- 
ance equation at the high electron temperature (Rees et al. 
1982). This makes everything rather complicated than the 
case of MPCS07. Conversely, once we know the properties 
of two-temperature plasma in MHD accretion flows which 
can reconstruct a radiation spectrum using a few relevant 
parameters, we can then constrain the physics of heating 
mechanism. 

In this study, we investigate radiative signatures of ra- 
diatively inefficient accretion flows in long-term 3-D global 
MHD simulations. Then, confronting with the observed 
spectra of Sgr A* in the flaring and quiescent states, we 
derive constraints for electron temperature that can concor- 
dant with the observations. We also discuss heating mecha- 
nism of electron in the magnetized accretion flows. In §2 we 
present our 3-D MHD model of accretion flows and describe 
method of radiation transfer calculation. We then present 
our results in §3. The final section is devoted to summary 
and discussion. 



2 NUMERICAL MODELS AND METHODS 
2.1 Setup of simulations 

A physical model we use in this study is based on 3-D re- 
sistive MHD calculations with pseudo-Newtonian potential 
(see Kato, Mineshige, Shibata 2004; Kato 2004 for more de- 
tails). We pick last 40 snapshots of our calculations from 
t = 30, 000 rs/c to 32,000rs/c with dt = SOr^/c where 
and c are the Schwartzschild radius and the speed of light, 
respectively. Note the flow has evolved in more than 600 ro- 
tation periods at the innermost stable circular orbits (ISCO) 
where the Keplerian rotation period is about 50rs/c. There- 
fore our MHD model is supposed to be in quasi-steady state. 
In our radiation transfer calculations, we use uniform meshes 
as (A'x, Ny, A^z) — (100, 100, 100) in Cartesian coordinates in 
a simulation box |a::,2/,2:| ^ 100 r-g. 



2.2 MHD model 

The quasi-steady MHD accretion discs have a hot, geomet- 
rically thick structure associated with sub-thermal magnetic 
held. Fig. [T] displays a spatial distribution of the density, the 
gas temperature, and the strength of magnetic held of our 
MHD model. The density is normalized by the initial maxi- 
mum density po and the strength of magnetic fleld is propor- 
tional to (p/po)^^^ for the same black hole mass A'Ibh- As we 
can notice, an MHD disc has non-axisymmetric structures 
(close to m = 1 where m is the azimuthal mode number) in 
density, gas temperature, and strength of magnetic field, si- 
multaneously. Moreover, filamentary structures of cold and 
dense gas can be seen in left and middle panels. Note that 
gas temperature is relatively high at funnel region along the 
z-axis because the centrifugal barriers prevent the penetra- 
tion of accreting material. In right panel, MHD turbulence 
induced by MRI and the differential rotation in the ffow 
generate strong magnetic field regions within approximately 
30r,. 



2.3 Electron temperature 

In our MHD model, it is assumed that the entire magnetic 
energy released in the diffusion region is thermalized in- 
stantly and therefore the production of non-thermal elec- 
trons accelerated by the magnetic reconnection is neglected 
for self-consistency. In the previous study OKM05, electron 
temperature is determined by the local thermal equilibrium 
between radiative cooling and electron heating via coulomb 
coupling in the cylindrical region. However, this assump- 
tion is invalid. Actually, in our preliminary radiation trans- 
fer calculation coupled with the energy balance equation, we 
found that the heating rate of coulomb coupling cannot af- 
ford the cooling rate of electrons in each computational cell. 
Therefore, the other heating mechanism of electrons, such as 
turbulent heating, must be taken into account for physical 
reasoning. Similar conclusions have been made by Sharma, 
Quataert, Hammett, & Stone (2007: hereafter SQHS07). For 
this reason, we introduce a new parameter, /op = Tc/Tp, 
the ratio of electron temperature and proton temperature so 
that electron temperature is determined with the gas tem- 
perature, Tgas, by using /op in this study. We assume that 
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/op is spatially uniform for simplicity. Thus electron temper- 
ature is obtained as: 

To = min ( f"" Tgas.mocVfcb ) (1) 

\ t + ,/op J 

where Tgas is the gas temperature which is derived by MHD 
simulation. We assume here that the electron temperature 
cannot exceed the temperature of rest mass energy of elec- 
tron due to pair annihilation. 




2.4 Parameters for radiation transfer calculation 

There are only three model parameters, po, /op, and Mbh, in 
order to perform radiative transfer calculation in our study. 
The mass of SMBH in the GC is fixed at Mbh = 3.6x lO'^M© 
(e.g., Schodel et al. 2002). The other model parameters 
we choose here are as follows: (a) po = 8 x lO^^^gcm""^ 
and fcp = 0.25, (b) po = 8 x lO^^^gcm"^ and fcp = 1, 
(c) po = 8 X 10~^®gcm-^ and /op = 0.25, and (d) po = 
8 X 10~^*gcm~^ and /op — 1. Note that these parameters 
are derived by fitting the observed broadband spectra. Be- 
cause the flow velocity is identical in all models, the larger 
density represent the more mass accretion rate A/. The rela- 
tion between po and M at ISCO can be described as follows: 

Mss2.5xl0"M ^ -|M0yr-\ (2) 

\^8 x lO-i^gcm-3 J '^^ ^ ^ 

This relation indicates that mass accretion rate of our mod- 
els is smaller than that for the quiescent X-ray emission 
measured with Chandra of M ~ 10~^ Mq yr~^ at the Bondi 
radius (Baganoff et al. 2003), but is consistent with that es- 
timated by the Faraday rotation in the millimeter band of 
M ~ lO"'^ - 10"* Mq yi-^ (Bower et al. 2003, 2005). 

In radiation transfer calculation, we treat synchrotron 
emission/absorption, (inverse-) Compton scattering, and 
bremsstrahlung emission/absorption of the thermal elec- 
trons. Non-thermal electrons produced by collisions of pro- 
tons via TT-decay (Mahadevan et al. 1998) are not taken into 
account, because we focus only on the radiative properties 
of thermal electrons (Loeb & Waxman 2007). 



Figure 2. Geometry of radiation fields, I in the coordinate space 
(x,y,z,6,<p) and coordinate screen {u,v) for the ray-tracing cal- 
culation. 

2.5 Radiative transfer calculation 

We solve the following radiation transfer equations with elec- 
tron scattering by using Monte-Carlo method: 

n-VI^^xAS.-Ir.) (3) 

where Iv(x,y,z,9,(l)) is the specific intensity at the posi- 
tion [x, y, z) in the direction {6, (f>) with the frequency v, 
whereas Xv{x,y, z, 9, (f>) = UcOv -f hii, is the extinction coef- 
ficient where no, (Ti, and is the electron number density, 
the scattering cross-section, and the absorption coefficient, 
respectively, and 

Sv = h Q> (p{v,n;u' ,n')a^iI,y{n')dQ,' (4) 

is the source function where Su, a^, and 'p{v, n; v'n') is the 
local emissivity, the scattering albedo, and the phase func- 
tion, respectively. The coordinate system is shown in Fig. 

El 

For generating photon packets, we randomly selected a 
position by using a local emissivity as follows: 

k— 1 Nmosh k 

where £iy(i) is the emissivity at the position index i and Ri 
indicates a random number distributed uniformly in the in- 
terval [0, 1]. In our study, all pseudo- random numbers i?j 
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are generated by using a Mersenne Twister method (Mat- 
sumoto & Nishimura 1998). The direction of generated pho- 
ton packets, n = (sin 6^ cos sin sin cos 6), is also de- 
termined by using random numbers R2 = {cos 6 + l)/2, 
Ri ~ (i>/2n. The frequency domain of photon packets is 
ranging from 10^ to lO'^^ Hz and is uniformally divided by 
100 bins in logarithmic scale. In this study, A'^p = 10^ pho- 
ton packets are generated in every frequency bin in order to 
acquire the statistically significant results. 

In order to evaluate an escaping probability of an emerg- 
ing photon packet, optical depth is computed by using a 
direct integration along the photon packet trajectory as fol- 
lows: 



(ncCTv + i^v) ds, 



(6) 



where / is a distance between the origin of the photon packet 
and the computational boundary along the trajectory. Here, 
we use the Klein-Nishina formula of scattering cross-section 
o"KN for (T^ (Rybicki & Lightman 1979) and the synchrotron, 
free-free, and bound-free self-absorption coefficient for Ki, de- 
scribed by Kirchoff 's law assuming local thermal equilibrium 
(LTE) at every meshes. 



+ £g+£,°' 



(7) 



where e^J , and are synchrotron, free-free, and bound- 
free emissivity, respectively (Pacholczyk 1970; Stepney & 
Guilbert 1983), and is the Planck function. Accordingly, 
an escaping probability of the generated photon packet is 
written as: 



w{V) = exp (-r^(Z)), 



(8) 



and the remaining photon packet 1 — 10(1) interacts with 
gas. Scattering position at the distance Is from the original 
position is determined by using a random number as follows: 



7?4 = [1 - w{k)] I [1 - w{l)\ 

and the scattering albedo is given by: 



(9) 



(10) 



where (kv) is the mean absorption coefficient along the pho- 
ton packet trajectory. To determine the phase function for 
the Compton scattering process, <p{y,n\v' ,n'\ we follow 
the method by Pozdnyakov et al. (1977). We repeat the 
same procedure for a scattered photon packet until either 
a photon is out of the computational box or wiV) < e where 
e = 10~^ in this study. Note that photons can neither pen- 
etrate nor be generated at the region within 2 Vs in above 
calculation. 

Finally, the computed radiation field is accumu- 
lated in the data array of {Ny^, Ny, N,., Ng , N^, N,j) = 
(100, 100, 100, 3,4, 5) in order to generate synthetic images 
in a given frequency band and viewing angle. On the other 
hand, the computed spectrum of escaping photons in the 
fixed viewing angle is stored in the different data array. 



3 RESULTS 

3.1 Mapping of escaping photons 

We investigate the spatial distribution of emergent radi- 
ation in order to explore the radiative nature of magne- 
tized accretion flows. Fig. |3] represents synthetic images 
of our models (a), (b), (c), and (d) at different frequency 
bands with the viewing angle of {0, <j))={n/3, 7r/4). In Fig. 
[31 upper half images correspond to the model with high- 
density plasma {M ~ 2.5 x lO~'^Af0 yr~^), whereas lower 
half images correspond to that with low-density plasma 
(M ^ 2.5 X 1O"*M0 yr"^). Each density model has two sub- 
categories; one is two-temperature plasma /cp — 0.25 and 
the other is one-temperature plasma fep = 1. A basic fea- 
ture of all images is that the core of emission region is larger 
than the size of the event horizon shadow, and is smaller 
than the radius of 50ra, except for the image of model (b) 
(namely high-density and one-temperature plasma) at mil- 



limeter band 10 



10"" Hz. 



Distinctive feature of non-axisymmetry is seen in all 
models at millimeter band. The non-axisymmetric structure 
in the magnetized accretion flows is also visible in Fig. [T] 
The apparent size of emission region in models (a) and (d) 
looks similar, whereas that in models (b) and (c) looks quite 
different. This is because the difference of density and tem- 
perature between model (a) and (d) compensate each other. 
On the other hand, non-axisymmetric feature cannot be seen 
in sub-millimeter band (10^^ — 10^^ Hz) except for model (b), 
but it has nearly spherical emission region around the grav- 
itational centre. Again, the apparent size of emission region 
in model (a) and (d) looks quite similar. At both millimeter 
and sub-millimeter bands, the core of emission region cor- 
respond to the region of the strong magnetic field (see Fig. 
[TJ;), and the peak of emission region is slightly shifted from 
the gravitational centre. 

All models represent the disc-like structure in X-ray 
band lO'^^ — 10^® Hz. The apparent size of emission region in 
high-density models [(a) and (b)], and in low-density mod- 
els [(c) and (d)] looks similar each other. This is because X- 
ray emission is produced by the Bremsstrahlung emission, 
which is sensitive to density, not to electron temperature. 
Interestingly, asymmetric structure can be seen only when 
electron temperature is equal to proton temperature [mod- 
els (b) and (d)]. This feature is most likely to be produced 
by Compton scattering of low-energy photons generated by 
synchrotron emission (known as synchrotron self-Compton 
process: SSC). 



3.2 Spectral energy distribution (SED) 

In addition to the spatial distributions of radiative fields, 
we calculate the SED of our MHD model. In Fig. |4l the 
resultant SED is shown from radio to gamma-ray bands 
and also the observed spectra of Sgr A* are superimposed. 
The overall spectrum consist of several radiative processes, 
that is, self-absorbed synchrotron for u < 10^^ Hz, syn- 
chrotron for lO^'^ < < 10^** Hz, synchrotron self-Compton 
for 10^'' 10^^ Hz, and thermal Bremsstrahlung for 
10"'^ Hz. As we can see, time variability is different at 
each band. In optically thick region below the critical fre- 
quency Uc ~ 10^^ Hz, there is a small time variation. When 
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Figure 3. Synthetic images of four models [(a), (b), (c), and (d) from top to bottom] in different frequencies (left: 10^*^ ^ < 10^^ 
[Hz]; centre: lO" < lO^^ [Hz]; right: lO^'' < W^^ [Hz]) at the viewing angle of (61, (/))=(7r/3, 7r/4). Colours indicate the power of 
escaping photon packets. A dotted circle with the radius 5rs indicates the canonical size of event horizon shadow. 
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the frequency becomes higher, close to the critical frequency, 
the time variation becomes larger. Note that the variability 
at the lowest frequency edge ~ l(f Hz depends on the num- 
ber of photon packets used in the calculation. Because of 
the strong self-absorption of synchrotron process, we can- 
not exclude the statistical errors in those frequency bands. 
Estimating from higher frequency regions, the variability is 
about two factors in magnitude. In optically thin region, 
on the other hand, time variation becomes prominent in the 
low-frequency band 10^'^ — 10^*^ Hz and it becomes negligible 
in high-energy band > 10^^ Hz. 

The resultant SED of model (a) can nicely fit the ob- 
served SED in flaring state for all frequency range (Fig.|3^). 
The difference between the emergent spectra and the total 
emissivity spectra below the critical frequency is induced by 
the synchrotron self-absorption, whereas that in X-ray band 
is caused by the effect of viewing angle. Unlike the previous 
study of ADAF and OKM05, low- frequency bump in the 
millimeter band 10^" — 10^^ Hz is reconstructed successfully 
without considering non-thermal electrons. This is a main 
product of introducing a parameter /cp, which can account 
for spatially structured electron temperature distribution. 

In order to explain two orders of magnitude difference 
in X-ray emission between flaring and quiescent states, the 
density need to be reduced by at least one order of magni- 
tude, such as model (b) (see Fig. |4}3). The resultant SED 
of model (b) can fit the X-ray spectra in quiescent state. 
However it under-predicts the radio to IR spectra for two- 
temperature plasma with /op = 0.25. Remaining option is 
to increase electron temperature so as to compensate the 
reduction of density, such as model (d) of one-temperature 
plasma with /ep = 1 (see Fig. |4ji). It turned out that the 
resultant SED of model (d) can successfully reconstruct the 
X-ray spectra of quiescent state as well as the radio and IR 
spectra. Although the underlying physics of electron heating 
is not clear, this is an outstanding result for the theory of 
hot accretion flows. 

It is interesting to note that the X-ray variability does 
not strongly depend on the model parameter. Actually, all 
models except model (b) represent that the X-ray variability 
is very weak. In model (b) , soft X-ray photons « 10^^ Hz are 
strongly contaminated by SSC photons because both density 
and electron temperature are increased so that scattering co- 
efficient becomes large. Moreover, it seems that most of vari- 
abilities in ultraviolet band are produced by the scattered 
photons and also they are correlated with IR variabilities 
generated by synchrotron emission in all models. Note that 
the amplitude of variabilities (different between the maxi- 
mum and the minimum power of escaping photons) in IR 
and UV bands does not change much in the range of our 
model parameter. Therefore, we suspect that the variabil- 
ity in the recent NIR observation (Eckart et al. 2005) are 
strongly affected by either the structure or the dynamics of 
the accretion flows. 



of the event horizon shadow, and is smaller than the radius 
of 100 Ts. We have found that a non-axisymmetric structure 
of m = 1 is associated with an elongated core emission in 
millimeter band. Remarkably, the peak location of the core 
emission in millimeter band ~ lO'^^ Hz is slightly shifted 
from the gravitational centre. This is consistent with the 
baseline-correlated flux density diagram of the recent VLBI 
observation at 230 GHz (Doeleman et al. 2008). 

We have also demonstrated for the first time that our 
3-D MHD model with different density (namely different 
mass accretion rate) can reconstruct the observed broad- 
band spectra including both the X-ray quiescent and flaring 
states, simultaneously. We have found that the X-ray flaring 
state corresponds to relatively high-mass accretion rate with 
a weak coupling between electron and proton temperature 
/cp = 0.25, whereas the X-ray quiescent state corresponds 
to relatively low-mass accretion rate with a strong coupling 
between electron and proton temperature /ep = 1. This is 
an opposite sense if one considers only the Coulomb cou- 
pling for electron heating. We will discuss this issue in the 
following. 

Heating mechanism of electron in hot accretion flows 
has been investigated by many groups (Bisnovatyi-Kogan & 
Lovelace 1997; Quataert 1998; Gruzinov 1998; Blackman 
1999; Quataert & Gruzinov 1999; Medvedev 2000). In 
ADAF/RIAF models, the turbulent viscosity primarily 
heats protons and then hot protons heats electrons via the 
Coulomb coupling between them. When the proton temper- 
ature becomes the virial temperature « lO^'^ K, the electron- 
proton Coulomb relaxation time becomes much larger than 
the dynamical time (Spitzer 1956; Stepney 1983). As a re- 
sult, the electron temperature decouple from the proton 
temperature. This is a well-known understanding of physics 
in hot accretion flows (Rees et al. 1982; Narayan et al. 1995). 
However, the spectra of hot accretion flows modeled by 
MHD simulations indicate that a coupling ratio /op is close 
to unity in both the flaring and quiescent state. Moreover, 
/op increases when the mass accretion rate decreases. This 
is an opposite sense in terms of plasma physics, because the 
electron-proton Coulomb relaxation time increases when the 
density decrease as ireiax oc Te^'^/po- Therefore, our results 
imply that the alternative heating mechanism of electrons 
is requisite in order to keep the electron temperature being 
close to the proton temperature. 

Recently, SQHS07 have found that signiflcant fraction 
of dissipative energy generated by turbulent viscosity can 
be directly received by electrons. This can naturally explain 
the discrepancy of /op between our results and ADAF/RIAF 
models, because a fraction of energy received by electron is 
assumed to be typically 5 ~ ma/m-p ~ 10"'' in ADAF mod- 
els (e.g., Narayan et al .1998). Although the detailed physics 
on the basis of viscous heating of electrons is in dispute, their 
results are worth to implement as subgrid physics of electron 
heating for 3-D global MHD simulations in the near future. 



4 SUMMARY AND DISCUSSION 

We have investigated three-dimensional radiative features 
of radiatively inefficient accretion fiows modeled by MHD 
simulations. The synthetic images of all models show that 
the core of emission region in Sgr A* is larger than the size 
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Figure 4. Spectral energy distributions for a different plasma model, (a) po 
10" 



8 X 10"-^^ [gem- 

'[gcm-3] and /ep = 1, (c) po = 8 X lO"!'' [g cm-^] and /ep = 0.25, and (d) po = 8 x 10-i'*[{ 



and /ep = 
m^"^] and fc 



0.25, 
, = 1 



(b) PO = 8 X 
A black solid 



line and a grey dashed line indicate the emergent spectra and total emissivity spectra, respectively. Black bars associated with solid and 
dashed lines indicate the time variation at a given frequency band. Blue crosses and arrows indicate data taken by radio and infrared 
(IR) observations (Narayan et al. 1998 and reference therein). Green bars indicate data by near infrared observation (Eckart ct al. 2006). 
X-ray observations of the flaring and the quiescent state by Chandra are depicted as red and orange "bow ties" , respectively (Baganoff ot 
al. 2001) whereas gamma-ray observation of the flaring state by INTEGRAL is shown in light-green "bow ties" (Belandgor ot al. 2005). 



diation transfer code have been developed on FIRST simu- 
lator at the center for computational sciences, University of 
Tsukuba. Radiation transfer computations were carried out 
on FIRST simulator and XT4 at the center for computa- 
tional astrophysics (CfCA), National Astronomical Obser- 
vatory in Japan (NAOJ). MHD computations were carried 
out on VPP5000 at the Astronomical Data Analysis Center 
of the National Astronomical Observatory, Japan (yyk27b). 
This work was supported in part by the FIRST project 
based on Grants-in-Aid for Specially Promoted Research 
by MEXT (16002003, MU), Grant-in-Aid for Scientific Re- 
search (S) by JSPS (20224002, MU), and Grant-in-Aid for 
Scientific Research of MEXT (20740115, KO). 
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